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Abstract. — We prove that the Navier-Stokes, the Euler and the Stokes equations admit a 
Lagrangian structure using the stochastic embedding of Lagrangian systems. These equations 
coincide with extremals of an explicit stochastic Lagrangian functional, i.e. they are stochastic 
Lagrangian systems in the sense of [5]. 
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1. Introduction 

The current paper aims at finding out a Lagrangian structure for some partial differential 
equations including the Euler, the Stokes and the Navier-Stokes equations. By a Lagrangian 
partial differential equation we mean that we can be find an explicit Lagrangian function and 



2 



JACKY CRESSON & SEBASTIEN DARSES 



an adapted variational calculus, for example the quantum, fractional or stochastic calculus of 
variation, such that the solutions of the PDE are critical points of the Lagrangian functional. 
This is the classical inverse problem for partial differential equations. 
The interest of such results is at least twofold: 

- First of all, we now have an intrinsic object, the Lagrangian, which control the dynamical 
behaviour of the PDE. 

- Secondly, this intrinsic structure can be used to derive more adapted numerical schemes 
for these equations. This will be discussed in a forthcoming paper. 

However, classical results in the study of incompressible fluids ensure us that there exists 
obstructions to obtain a variational principle from which the behaviour of the fluid derives. 
An analogous result exists for the derivation of a Lagrangian structure for non conservative 
mechanical systems (see Bateman [1]). However, these obstructions results are based on the 
fundamental assumption that we use the classical differential calculus. As a consequence, 
these obstructions may disappear if we use a different "differential" calculus by extending the 
underlying class of objects to stochastic processes for instance. 

Stochastic interpretations of the Navier-Stokes equations already exit in a huge literature. 
For instance, C. peskin gives in |14) a random- walk interpretation based on a model where 
the force exerted by a molecule of the fluid is considered. In [2] B. Busnello turns the Navier- 
Stokes equations into a system of functional integrals in the trajectory space of a suitable 
diffusion. In [3j, the authors provide a stochastic representation where the drift term in a 
stochastic differential equation is implicitly computed as the expected value of an expression 
involving the flow it drives. This is a step towards a generalization for the Feynman-Kac 
theory for parabolic equations. 

In the current paper, we use the stochastic embedding of Lagrangian systems developped 
in [41, [6] to relate the Navier-Stokes, the Euler and the Stokes equations to critical points 
of a stochastic Lagrangian. An embedding theory is the data of a functional space and an 
extension of the classical derivative on this set. This allows us to define natural extension 
of differential operators and differential equations. When dealing with Lagrangian systems 
an embedding allows us to keep track of the form of the equation but also of the underlying 
variational principle. Using embedding, we then obtain a strong relation between a classical 
Lagrangian dynamics and an embedded one, which will be the PDE in our case. 

The body of our paper is arranged as follows: We first remind definitions and results about 
the Nelson derivatives for stochastic processes and the stochastic embedding of Lagrangian 
systems developped in [6]. We then prove that the Navier-Stokes, the Euler and the Stokes 
equations coincide with extremals of stochastic Lagrangian systems. 



THE STOKES, NAVIER-STOKES AND EULER EQUATIONS 



3 



2. Notations and Reminder about Nelson stochastic derivatives and Time 

reversal 

2.1. Notations. — Let T > Q, v > and d G N*. Let L be the set of all measurable 
functions / : [0,r] x M*^ — M'' satisfying the following hypothesis: There exists K > such 
that for all x,yeR'^: supj/(i,x) - f{t,y)\ <K\x-y\et sup^ < K{l + \x\). 

We are given a probability space {i},A,¥) on which a family (W^'-''''^^)(6,cr)eLxL of Brownian 
motions indexed by L x L is defined. If 6, cr G L, we denote by the natural filtration 

associated to W^^'"\ Let V be the filtration generated by the filtrations where {b,a) € 

L X L, and we set: For t G [0,T], 

(&,CT)eLxL 

Let F([0,T] X M'^) be the space of measurable functions defined on [0,T] x M"^ and let 
F([0, T] X $7) be the space of measurable stochastic processes defined on [0, T] x $7. 

Let us define the involution ^ : F([0, T] xM'^) F([0, T] xM'^) such that for aU t G [0, T] and 
X G M'', {(fm){t,x) = —u{T — t,x). We also define the time-reversal involution on stochastic 
processes: r : F{[0,T] x Q) ^ F{[0,T] x O), r{X)t{u) = XT-t{uj). 

It is convenient to use the bar symbol to denote these two involutions. We now agree to 
denote deterministic functions by small letters and stochastic processes by capital letters. So 
there will not be any confusion when using the bar symbol: u := (f>u and X := r{X). 

The space is endowed with its canonical scalar product (•,•). Let | • | be the induced 
norm. If / : [0,r] X M'^ ^ R is a smooth function, we set djf = We denote by 

V/ = {dif)i the gradient of / and by A/ = '^jdjf its Laplacian. For a smooth map 
$ : [0,T] X R"^, we denote by its j*^-component, by (dx^) its differential which we 

represent into the canonical basis of M'^: (dx^) = {dj^^)ij, and by V • $ = Y2j dj^-' its diver- 
gence. By convention, we denote by the vector (A^^)j. The notation ■ V)<I> denotes 
the parallel derivative of <I> along itself, whose coordinates are: ((<!> • V)<I>)* = ^Wj^^. 
The image of a vector u G M'^ under a linear map M is simply denoted by Mu, for instance 
(dx^)u. A map a : [0,T] x M*^ M'^ (g) M'^ is viewed as a d x d matrix whose columns are 
denoted by a^. Finally, we denote by V • a the vector (V ■ ak)k- 

By the Navier-Stokes equation we mean the following equation: 



(1) 



dfU + {u ■ V)n = vAu — Vp, 
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where u : [0,T] x M*^ ^ M'^ is the velocity field and p : [0,T] x M'^ ^ M is the pressure. 
By the Navier-Stokes equations we mean the following system: 

(2) dfU + (n • V)u = i^Au - Vp, V ■u = 0. 

When talking about solutions of ([T]) (resp. ([2|)), we will always refer to a pair (n,p) of strong 
solutions of ([T|) (resp. ([2])). We distinguish the equation ([T|) from the system ^ because 
in the current paper we are interested in the structure of ([T]). The system 1^ turns out to 
have the same Lagrangian structure but in that case we have to embed this Lagrangian into 
a restricted class of diffusions processes, namely those with a divergence free drift. 
We agree to use the same convention for the Euler and the Stokes equations. 

2.2. Time reversal and Nelson derivatives. — 

Definition 1. — We denote by the space of all diffusions X satisfying the following 
conditions: 

(i) X is a solution on [0,T] of the SDE: dXt = b{t, Xt)dt + a{t, Xt)dwf '''^ Xq = X^ 
where X^ e L^(0) and (6, cr) G L x L, 

(ii) For all t £ (0, T), Xt admits a density pt{ ) , 

(iii) Setting a*-^ = {aa*y^ , for all i G {1, • • • , n}, to > 0, 

\dj{a^\t,x)pt{x)) \ dxdt < +oo, 



to -'K'* 



(iv) For all i,j, t, 



d,{a'\t,.)pt{-)) 



G L. 



A diffusion verifying (i) will be called a (6, cj)-diffusion and we denote it by X^^''^^ 

We denote by the closure of Vect(A-^) in L^{^} x [0,T]) endowed with the usual norm 

We know that the reversed process of any element A^ is still a Brownian diffusion driven 
by a Brownian motion W^^''^'^ {cf |12] ). We denote by V^^'"^^ the natural filtration associated 
to Let V be the filtration defined by: For ah t G [0,T], 



(fe,o-)eLxL 



(b,a) 
t 



We finally consider the filtration such that J^t = Vt-i for all t G [0,T]. 



Proposition 1 (Definition of Nelson Stochastic derivatives) 
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Let X = X^^" G with a'^ = {aa*)'^ and = (a^^ • • • ,a'^^). For almost all t G (0,r), 
the Nelson stochastic derivatives exist in L'^{Q) : 

Xt+h — Xt 



(3) DXt : 

(4) D,Xt : 



lim E 
lim E 

h^0+ 



h 

Xt - Xt- 



t-h 



h 



Vt 



b{t,Xt) 

b{t,Xt) -^^Y.dj{a^{t,Xt)pt{Xt)). 



Therefore, for X^''^ G A^, there exists a measurable function 6* such that D^:Xt = h^{t,Xt). 
We call it the left velocity field of X. It turns out to be an important objet for the sequel. 
Also, this field is related to the drift of the time reversed process X through the following 
identity: 



(5) 



DX = h. 



We denote by = {X £ A^; DX, D^X G A^} and A^ the closure of Vect(A2) in L^{9. x 
[0, r]). We then define A^' in an obvious way. 

We denote by T)^ the stochastic derivative introduced in ([4j Lemme 1.2) and defined by 



(6) 



D + D^ D-D^ 
V„ = : + ^ /i G {0,±l,±i}. 



2 " 2 

We can extend V by C-linearity to complex processes A^ := Aj © iA^. 

Theorem 1. — Let X^^-'^) g A^ a = aa* and f G C^'^{I x M"^) such that dtf, Vf and dijf 
are bounded. We get: 



(7) 



Moreover, we can generalize for the operator the "product rule" given by Nelson in |13j 
p. 80, which is a fundamental tool to develop a stochastic calculus of variation: 

Lemma 1. — Let X,Y £ A}.. Then E[V^,Xt -Yt + Xf V^^Yt] = jiE[Xt ■ Yt]. 

The proof is a an immediate consequence of the form of the operator T) and to the fact 
that A"^ is a subspace of the class S{J-,Q) ([17J p. 226) for which W. Zheng and P-A. Meyer 
show the product rule of Nelson (c/ [l7j Th. 1.2 p. 227). 



3. Reminder about stochastic embedding of Lagrangian systems 

We first define the stochastic analogue of the classical functional. 
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Definition 2. — Let L be an admissible Lagrangian function. Set 

f-T 

\L{Xt,V^Xt)\dt <oo 





The functional associated to L is defined by 



(8) 



F : 



X 



E 



LiXt,V^Xt)dt 



In what follows, we need a special notion which we will call L-adaptation: 

Definition 3. — Let L be an admissible Lagrangian function. A process X £ is said to 
be L- adapted if: 

(i) X G E; 

(ii) For all t G /, d^L{Xt,V^Xt) G L^{n); 

(iii) d,L{Xt,V^Xt) G 

The set of all L— adapted processes will be denoted by C 



We introduce the following terminology: 

Definition 4- — Let T be a subspace of and X G A^. A T -variation of X is a stochastic 
process of the form X + Z , where Z £T. Moreover set 

r= = {z G r,VA: g + x g h}. 

We now define a notion of difjerentiable functional. Let F be a subspace of A^. 

Definition 5. — Let L be an admissible Lagrangian function and F the associated func- 
tional. The functional F is called T-differentiable at X £ C if for all Z £T-= 

(9) F{X + Z) - F{X) = dF{X, Z) + R{X, Z) , 

where dF{X, Z) is a linear functional of Z £T-= and R{X, Z) = o{\\ Z \\). 

A T-critical process for the functional F is a stochastic process X £ 3 D C such that 
dF{X, Z) =0 for all Z G Fh such that Z{a) = Z{b) = 0. 

The main result of [6] is the following analogue of the least-action principle for Lagrangian 
mechanics. 

Theorem 2 (Global Least action principle). — Let L be an admissible lagrangian with 
all second derivatives bounded. A necessary and sufficient condition for a process X £ CCiA^ 
to be a A^ -critical process of the associated functional F is that it satisfies 



(10) 



dL 

dx 



{Xt,V^Xt)-V. 



dL 

dv 



iXt,V^Xt) 



0. 
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We call this equation the Global Stochastic Euler- Lagrange equation (GSEL). 



The main drawback of Theorem [2] is that equation (jlOp is not coherent (see [6], §.6. 2) i.e. 
that it does not coincide with the direct stochastic embedding of the classical Euler-Lagrange 
equation of the form 



(11) 

except when /x = 0. 



dL_ 
dx 



{Xt,V^Xt)-V, 



dL 

dv 



{Xt,V^Xt 



0. 



In order to obtain a coherent embedding without imposing ^ = 0, we must restrict the set 
of variations. Let us introduce the space of Nelson differentiable processes: 

(12) = {X e h^.DX = D^X}. 

Using N^-variations we have been able to prove the following result [6] : 

Proposition 2. — Let L he an admissible lagrangian with all second derivatives bounded. A 
solution of the equation 



(13) 



dL_ 

dx 



{XuV^Xt)-V, 



dL 

dv 



iXt,'D^Xt) 



0, 



called the Stochastic Euler-Lagrange Equation (SEL), is aJ^"^ -critical process for the functional 
F associated to L. 



We have not been able to prove the converse of this lemma for ]M"^-variations. 



4. Stochastic embedding and the three equations of fluid mechanics 

We consider the Lagrangian function 



(14) L{t,x,v) = — - p, 

where p is a function depending on x and t. The classical Euler-Lagrange equation yields 

(15) ^ = -Vp. 

In the sequel, we investigate the stochastic embedding of this lagrangian over diffusion 
processes of specific forms. Precisely, the Navier-Stokes equation, the Euler equation and the 
Stokes equation turn out to be respectively related to the following subspaces of A^: 

1. The class composed of all the diffusions X*^"''^^ such that a is constant, the drift u is 
C^, bounded with all its second derivatives bounded, and Vlogpt has bounded second 
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order derivatives. We can prove that Kl C thanks to Prop. 2 p.394 in [7]. Recall 
that a (n, cr)-diffusion with a constant diffusion coefficient a is of the form: 

(16) Xt = XQ+ I u{s,Xs)ds + aWt. 

Jo 

2. The class Ag composed of all the smooth processes of the form: 

(17) Xt = Xo+ [ u{s,Xs)ds, 



Jo 

where u G C2(M+,]R'^) and Xq £ L^{n). 
3. The class A^ composed of all the processes of the form: 

(18) Xt = Xo+ [ u{s,aBs)ds, 

Jo 

where u G C2(M+,M'^), a £ R, Xq e L^{n) and B is a Brownian bridge pinned to be 
at T. 

5. The incompressible Navier-Stokes equation 

Recall that the incompressible Navier-Stokes equation can be written as 

(19) dtu + (n • V)^ = uAu - Vp. 

Regarding stochastic least action principles, we can consider two cases: The first one uses 
a full set of variation and the second one uses N^-variations. 

5.1. Full variations. — Applying Th. [2l a (u, (T)-diffusion X which is a critical point of 
the natural lagrangian (jl4p for full variations, satisfies 

(20) (V^Xt) = -Vp. 

When fi = 1, i.e. D = D, we obtain D^:{DXt) = D^u{t,Xt) = —Vp and finally: 

2 

(21) dtu + u^ -Vu- Y^u = -Vp. 

When n = —1, i.e. V = D^, we obtain D{D^Xt) = Du^{t, Xt) = —Vp and finally: 

(22) + u ■ Vu* + — An* = -Vp. 

We notice that the field u does not satisfy the Navier-Stokes equations when using full 
variations in the case ^ = 1 since the constant in front of Au is positive. When /i = — 1, 
we can obtain the Navier-Stokes equation if and only if u^: = u, i.e. if DXt = D^Xt which 
corresponds to Nelson differentiable processes. These processes can be described through a 
given partial differential equation (see [6j). In particular, a diffusion belonging to A^ with a 



THE STOKES, NAVIER-STOKES AND EULER EQUATIONS 



9 



constant diffusion coefficient cannot be a Nelson difFerentiable process. As a consequence, we 
have no hope to recover the Navier-Stokes equation into this setting. 

5.2. 3sf^-variations. — One of the main results of ^ states that the stochastic Euler- 
Lagrange equation reads in that case: 

(23) {V^Xt) = -Vp. 

When ^ = 1, i.e. V = D, we obtain D'^Xt = Du{t, Xt) = —Vp and finally: 

2 

(24) dtu + u ■ Vu + yAu = -Vp. 

We notice that up to the sign of the constant term we obtain the form of the Navier-Stokes 
equation. As cr^ is always positive we have a apparent obstruction to obtain the Navier-Stokes 
equation using the D-embedding. However this obstruction can be removed considering drifts 
of the form u. 

When fj, = —1, i.e. P = D^,, we can obtain the Navier Stokes equation from a {u,a)- 
diffusion and it will be satisfied by n*. 

These results are the content of the following theorem: 



Theorem 3. — Let u > and ^("'^^i.) G A^. 

If the left velocity field ti^,(t, x) of X satisfies the Navier-Stokes equation 
(25) dtu^, + n* • Vii* — i/An* = —Vp, 



then the stochastic process X^^'^^^ is a -critical process of the stochastic functional 



(26) X 



T 



L{Xt,D,Xt)dt 







„2 



where L is the natural lagrangian L{x, v) = ^ — p. 

Moreover, if u verifies the Navier-Stokes equation then the stochastic process X^"'^^) is a 
-critical process of the stochastic functional 



(27) X^E 



L{Xt,DXt)dt 







,,2 



where L is the natural lagrangian L{x, v) = ^ -\- p. 
Proof — Let X^"'"^) G A^ and let be its left velocity field D^X. Then we have: 
(28) Z)*(D,A:f"''"^) = (^dtu^ + • Vu^ - yAn,^ (t,x("''^)). 

Therefore iiu^{t,x) satisfies the Navier-Stokes equation (125]) . then: 

D*(i?,X,("'^)) = -Vp(t,x("''^)), 



10 



JACKY CRESSON & SEBASTIEN DARSES 



SO X^^'^^^ is a !N""^ critical process of the stochastic functional 



(29) 



X 



L{Xt,D,Xt)dt 



Now, assume that u satisfies the Navier-Stokes equation and let us consider the diffusion 
X := Then 

(30) D'^Xt = {dtu + u-Vu + uAu){t,Xt) = {dtu + u-Vu- iyAu){T - t,Xt) 

(31) = -Vp{T-t,Xt)=Vp{t,Xt), 

which yields the last statement of the Theorem. 

□ 

Remark 1. — Due to the fact that Theorem [2] does not give an equivalence between J^^- 
critical process and the stochastic Euler-Lagrange equation, we are not able to prove that the 
Navier-Stokes can be obtain by a coherent embedding procedure, i.e. that the following 
diagramm commutes 

Emb(D,) 



(32) 



L{x{t),x{t)) 
LAP 



d dL 



dL 



^L{Xt,D,Xt) 

Stochastic LAP 



where z{t) = {x(t),x{t)), Z{t) = (X(,I?*X(), and LAP stands for "Least Action Principle". 



6. About the incompressible Euler equation 

The incompressible Euler equation 
(33) dtu + {u ■ V)n = -Vp 

is of course obtained as a limit of the Navier-Stokes equation when the viscosity i/ goes to 
zero. However, in our case this has a strong consequence on the underlying embedding from 
the functional point of view. Indeed, when v goes to zero, the set Sq reduces to deterministic 
stochastic processes which are sufficiently regular. On the contrary the basic property of 
diffusion processes is that each trajectory is a.s. Holderian with exponent 1/2 — e for e > 0. 
First of all, let us remark that D = ~ Tf on 5o. 



Corollary 1. — A process Xt G Sq of the form 

(34) dXt = u{t, Xt)dt 
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is a 3sr^ or -critical process of the stochastic functional 

r-T 



(35) 



X ^ E 



L{Xt,DXt)dt 



where L{x, v) = — — p, if and only if the velocity field u{t, x) satisfies the Euler equation 



(36) 



dtu + u ■ Vu = — Vp. 



One must be carrefull with this result. Even if the basic underlying set of stochastic 
processes corresponds to deterministic processes, the stochastic calculus of variations embed 
these deterministic processes into arbitrary diffusion processes. 



7. About the incompressible Stokes equation 

Through processes in the class As, we can find the following stochastic least action formu- 
lation of the Stokes equation: 

Theorem 4- — Let X ^ Kb he a process of the form: 

(37) dXt = u{t, y^Bt)dt. 
The velocity field u satisfies the Stokes equation 

(38) dtu - uAu = -Vp, 

if and only if the stochastic process Xt is a critical process of the stochastic functional 



(39) 



where L{x, v) = — — p. 



X^E 



r L{V2i)Bt,DXt)dt 
Jo 



Proof — The A^ least action principle reads: 

D^DXt = -Vp{V2^Bt). 

But 

(40) D^DXt = (dfu + \p2^ {d^u)D^Bt - i^Au^ {t, v^Bt). 

Since D^B = and the probability density of the bridge is everywhere positive, we can deduce 
the desired equivalence. 

□ 
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8. The temperature equations associated to the various fluid equations 

8.1. Temperature for the Navier-Stokes equation. — As mentioned in [15 j. one has 
to associate to the Navier Stokes equations an equation for the temperature field 9. This 
leads to the so-called Navier-Stokes-Fourier system (see e.g. |10]): 

(41) dtb + b-Vb- i^Ab = -Vp, V • 6 = 0, 

(42) dtO + b-Ve = kAO. 

The coefficient k is a priori different from the viscosity coefficient v. 

We have V • (pu) = u ■ Vp + pV • u. So the incompressibility condition V • 6 = implies 
that the temperature equation (|42p is actually the Fokker Planck equation: 

(43) dtp = -V • (pb) + kAp, 

satisfied by the probability density p of a diffusion x'^^'^^h 

Therefore the Navier Stokes system together with the temperature equation can be de- 
scribed by the two following diffusions: 

• X^'^'^^\ which is the critical point of the natural Lagrangian L{x,v) = '\ — p, and 

• x^"*'^^) whose density satisfies the associated temperature equation. 



8.2. Temperature for the Stokes equation. — The Stokes equation 

(44) dtu - vAu = - Vp, 

differs from the Navier-Stokes and Euler equations in the sense that it cannot be obtained as 
any limit of the previous ones. Moreover, the related equation on the temperature 9 does not 
surprisingly depends on the velocity field u: 

(45) dt9 = kA9. 

Based on Theorem [3l this specific structure might suggest to uncouple the underlying 
process of the velocity field and the velocity field itself. Based on this remark, we came up to 
Theorem HI Therefore the Navier Stokes system together with the temperature equation can 
be described by the two following processes: 

• Xt = Xq + Jq u(t,\/2uBt)dt, which is the critical point of the natural Lagrangian 

2 

-L(x, v) = ^ — p, and 

• \/2k B which is a Brownian motion and whose density satisfies the associated heat 
equation, that is the temperature equation ([I5]) . 
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8.3. Temperature for the Euler equation. — Corollary [T] naturally leads to the tem- 
perature equation for the Euler equation, thanks to elementary results about the transport 
of measure through an ordinary differential equation. 
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